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Abstract 

We consider the possibility to produce a bouncing universe in the framework of 
scalar-tensor gravity models in which the scalar field potential may be negative, and 
even unbounded from below. We find a set of viable solutions with nonzero measure 
in the space of initial conditions passing a bounce, even in the presence of a radiation 
component, and approaching a constant gravitational coupling afterwards. Hence we 
have a model with a minimal modification of gravity in order to produce a bounce in 
the early universe with gravity tending dynamically to general relativity (GR) after 
the bounce. 
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1 Introduction 


Since the realization that the strong curvature singularity arising in the past of our Universe 
in Friedmann-Lemaitre-Robertson-Walker (FLRW) models may not be avoided in generic 
solutions of the Einstein gravity with matter in the form of an ideal fluid (only its type 
changes from an isotropic one to the generic BKL vacuum singularity), there were many 
attempts to construct solutions avoiding this singularity by having at least one bounce 
at some high value of curvature either using more complicated held-theoretical models of 
matter or by modifying gravity. The simplest example to obtain one bounce in a set of 
solutions having a nonzero measure in the space of initial conditions is presented by a 
massive scalar held in a closed (i.e. positively spatially curved) FLRW universe p. This 
does not remove the curvature singularity in a generic solution, even in the class of FLRW 
models, but shifts it to the past, beyond any given hnite number of bounces of the scale 
factor. Non-singular FLRW solutions with an inhnite number of bounces exist, too, but 
they are degenerate: they have a zero measure in the space of initial conditions [2], see also 
[3]. This occurs because a minimally coupled scalar held with a non-negative potential can 
violate the strong energy condition, but not the weak and null ones0 However, having a 
bounce due to a positive spatial curvature requires severe hne tuning of initial conditions 
before the contraction stage in the past, see e.g. the calculations of the measure of those 
ones leading to a bounce in the above mentioned cosmological model in P and |6]. That 
is why it is much more interesting to obtain a non-degenerate bounce in FLRW models in 
the absence of spatial curvature. Previously known examples of such FLRW models are 
based on such radical modihcations of general relativity (GR) as theories with scalar P |H] 
or tensor ghosts, loop quantum gravity (see e.g. 0 ) or gravity described by an ehectively 
non-local Lagrangian (see e.g.|10]. HU, and [12] for a recent review). By contrast, in our 
paper we would like to restrict ourselves to the well-known and very modest modihcation 
of GR - scalar-tensor gravity without ghosts. On the other hand, we want to abandon 
the common assumption that the scalar held potential V(</>) in this gravity theory is non¬ 
negative, we will even consider the case when it is unbounded from below. In Section 2, 
the bouncing model is presented and general equations and results are given. In Section 
3, the existence of physically viable bouncing solutions is shown and their behaviour is 
studied in details. Section 4 contains our conclusions. 

^In the case of open (negatively spatially curved) FLRW models, a bouncing solution was found in [3] 
for a non-minimally coupled scalar field with a quartic potential (so this model represents a particular kind 
of scalar-tensor gravity models). However, later it was shown in [5] that this regular solution is degener¬ 
ate, too: it is unstable against arbitrarily small anisotropic and/or inhomogeneous metric perturbations 
preventing transition to the regime with a tensor ghost (repulsive gravity). 
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2 A bouncing model 


Let us consider a universe where gravity is described by scalar-tensor theory. The La- 
grangian density in the Jordan frame of the gravitational sector is given by 

L = i [F(<h)i? - Z(<h) - 2t/(<h)] . (1) 

Everywhere below we will use the freedom to take Z = 1 or Z = —1, corresponding 
physically to ojbd > 0 or ubd < 0 {ojbd = ZF (^) ^). For ojbd < 0, the theory is 
ghost-free provided —| < ubd < 0. We consider further spatially flat FLRW universes 
with metric ds^ = —dF -|- a?{t)dy? yielding the following modihed Friedmann equations 

- - 3iLF + t/ = 0 , (2) 

2FH + Z^‘^ + F-HF = Q. (3) 


with H = ^. Here and below a dot, resp. a prime, stands for the derivative with respect 
to t, resp. to <h. Equations (I2]),([3]) contain the equation of motion of <h 


Z (-$ - 3H^) + 3F'{H + 2H^) - U' = 0 . (4) 

We start our analysis with Z = 1 and we specialize to the following model 


F = + (5) 

0 

U = (6) 

where k~’^ > 0, A > 0 and c are constant parameters, only c being dimensionless. Sur¬ 
prisingly, while we arrived at this model by requiring that a combination of ([2]), ([3]) and 
dH), without derivation with respect to t, be zero, this model has been considered in the 
past in completely different contexts (see e.g. [13]). Thus (jl]) corresponds to a massless 
scalar held conformally coupled to Einstein gravity with an interaction potential — 0 $“^. 
Therefore, the energy-momentum tensor of the held $ is the sum of “dark” radiation with 
a negative energy density and a cosmological constant. 

As a result, the system ([2]), ([3]) reduces to 


1 

2 


3H‘^ 


k + K^ 


-4 



(7) 

( 8 ) 


with an ehective positive cosmological constant A and the bare gravitational constant 
SttG = = 0) = while y = a<h and g = f dt/a{t). The crucial point is that the 
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constant A - the energy density of the held x in Minkowski space-time - can be negative 
here. 

For y4 < 0, we have in front of ns an explicit illnstration of the possibility for scalar- 
tensor models to accommodate an effective dark energy component of the phantom type 
{Weff < -1) [I1],[I5] B Onr system can be completely integrated. As we are interested 
in bonncing solntions, we must take A < 0 and from (0), we have c > 0 so that U is 
necessarily an inverted potential, unbounded from below. In spite of looking unphysical 
at hrst sight, scalar helds with such an interaction have been often considered both in 
quantum held theory and cosmology, see e.g. m 

Integrating ([7]), a bouncing solution is obtained 


a = 

H = 



(9) 

( 10 ) 


where Oq = ^ ^ ^^e value of a at the bounce located at t = 0 with a trivial 

redehnition of t. It satishes H > 0 and has a constant Ricci scalar R = 6{H + 2H‘^) = 4A. 
Integrating (|8]) the analytical expression for x is obtained in terms of the Jacobi elliptic 
function dn(M|m) [IH] as follows 


xiv) 


^ 1 

c / dn(M|2) 


( 11 ) 


with 

where dn“^(M|m) is the inverse function of dn(M|m) and <Fo = = 0). The conformal 

time ri{t), chosen here such that //(O) =0, is given by 



where F{x\m) is the elliptic integral of the hrst kind [18]. This function is odd and then 
r]{t) is real. It is possible to show that 

1 T^fl) 

h(oo) = — 

phenomenological dark energy model in GR based on ([T]) was considered in m, leading to m and 

cni). 
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So we get finally 



We note that nnder the transformation Z —)■ — 1, F ^ —F ,U —U, the same solution 
flThll is obviously obtained and that ([7]) and (E]) remain valid as well because (I2]),(|3]) and 
(P|) are invariant under this transformation. 

Another attractive feature of this model is the possibility to have a bounce in the 
presence of a true radiation component resulting in a trivial change A —)■ A' in ([7]). This 
can model a bounce in the primordial universe. The bounce will survive provided the sum 
of “dark” and true radiation remains negative or 

A + Pradfi 0-0 = A' = -^ < 0 , (16) 

which shows that ^ = 4? is invariant, the scale factor a'n at the bounce decreases in the 

“o “o ’ ^ 

presence of radiation and also when A increases towards zero. Let us hnally note that the 
bounce disappears for A' > 0. 

3 Study of the bouncing solution 

In spite of the fact that we have obtained an explicit closed form expression for <h(f) we 
will study the behaviour of the bouncing solution by applying the qualitative theory of 
differential equations. This will provide us with a more intuitive and direct approach and 
will enable us to obtain results about the qualitative behaviour of our system otherwise 
difficult to hnd from ffTSj) . Solving for the wave equation recast into the form 

$ + 3H^ + 4c$(I'2 - $2) = 0 , (17) 

where $ = (^)T and evaluating (|2]) at the bounce 

li>S - c(4S - = 0 , (18) 

where <ho,mm = (^)^ = system is automatically solved. 

The inequality <1>2 < Qk~‘^ is readily obtained from the requirement F > 0 while $ = 0 
is impossible from ([H]). Choosing $o > 0, the following inequality is further obtained 

$ < ^0,min ^ *ho < ^max y (19) 
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Figure 1: The functions F (full) and U (dashed) are displayed for the model parameters 
= 20, c = 3 and A = 3. The two curves coincide at <F = 0 because we have actually 
plotted AF. While U vanishes at <F = <Fo,mm = F does so at <F = ^max = For 

these parameters we have d’o.mm = 0.4728, <Fo,cr = 0.4936860, ^max = 0.5477. For $ > 0, 
the range of possible values at the bounce lie in the interval <Fo,mm < 4* < ^max but only 
values in the interval ‘ho.cr < *h < *Fmax and with a negative slope will yield physically 
viable solutions. 


with 6 = ^max- Indeed, from (Ej), the condition 17(4>o) = Uq < 0 implies immediately 

*I’o > ^o,min, while <h < ^raax) Valid at all times, follows from the physical condition F > 0. 
A non vanishing interval allowed for $o requires <Fo,mm < ^max or 


A 


1 , 


( 20 ) 


implying also <F < From fflSj) . with <Fo < 0 (the justihcation of this choice will be 

clear below), we obtain 


4„ = -AcW-4S,».«). (21) 

with $0 decreasing when <Fo is increasing. 

We study the dynamics for f > 0. For f < 0 the dynamics follows immediately due to 
the invariance of the system under reflection t —t. We will come back to this point at 
the end of our analysis. We will analyse now the possible existence of local maxima and 
minima for ‘h(t). From (|8]) we readily obtain for $ = 0 


$2 


1^2 + - 4cAa-4 


( 22 ) 


The choice of the positive sign before the square root is compulsory since $2 must be 
positive. It is straightforward to show that fl20|) implies $ < ^max for $ given by ([22]). 
Hence the inequality <h < ^max follows when <F = 0. We have further from (Ej) when <F = 0 
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Figure 2: The function <F(t) is shown starting at t = 0 at the bounce for the same model 
parameters as in Figure 1. The initial value $0 is slightly below the critical value <Fo,cr = 
0.4936860... so this solution will diverge in a hnite time. 

that —4><F + QFH = 0, therefore $ > 0. So we have shown the impossibility for <F to have 
a maximum. 

On the other hand if <F has a hnite limit for t —)■ 00 , we must have <Foo = 0 and <Foo = 0, 
therefore flTTD gives 4c<Foo(®^ — = 0- Hence we have either $00 = 0 or $00 = *h. 

These results imply in particular that had we chosen $0 > 0, would eventually enter the 
unphysical domain <F > ^max, so this choice is not allowed. 

If the solution has a minimum, it is unique as <F does not have any maximum. This 
minimum must be larger that $ as a direct consequence of fflTj) . Hence after the minimum, 
<F > 0 forever since there is no maximum and then <F goes to inhnity - therefore crossing 
the value <hmax ^ because the only possible hnite limits $00 are lower than the minimum 
and cannot be reached. We will now show that in this case $ tends to inhnity in a hnite 
time too- When <F —)■ cxo, keeping the leading terms in ([2]) we easily hnd <F ~ and 

after integration 

$-pi-. (23) 

— y/^t + k 

with A; > 0 and too = As a consequence F becomes zero in a hnite time too, rendering 
this solution unphysical. When $0 = *Fo,min {Uq = 0), the bounce is a minimum from 
eqs. flTTl) and then <F goes to inhnity in a hnite time (see Figure [2]). 

Let us consider now the case where there is no minimum, i.e. <F < 0 for t > 0. From 
flTT)) we get <h = —3H^ — 4c<F(|)^ — <h^) > 0 as long as <h > I*. However when <F(t) < !>, 
the sign of <h(t) can change. Two subcases will therefore arise: 

1) <F(t) is always greater than $ and, since <h(t) < 0, tends asymptotically to = <F 
with positive concavity (‘h(t) > 0). 

2) <h(t) crosses <F and, since <F(t) < 0, tends asymptotically to 0. Below $, <F(f) can 
however change sign. 
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In order to differentiate these two possibilities, we study the system in the neighborhood 
of t = oo where ffTTj) can be recast as an autonomous planar system, viz. 




and A 2 = — 2 W hence it is a saddle point. So the solution $ —)■ d) is “unstable” and 


only with an exact initial value $0 can one reach the hxed point ($ 005 *^ 00 ) = (*h,0). In 


contrast, once we hnd some <I>o value which yields a solution tending to 0, we can vary $0 
in a neighbourhood of this value and still go to the hxed point ( 0 , 0 ). 

By continuity, if we take $0 slightly larger than *h(t) will still go to inhnity in 

a hnite time. Further increasing <Fo) two possibilities arise: 

a) For <Fo,mm < "ho < ‘ho.cr, where <I)o,cr is some critical value, ^(t) diverges in a hnite 
time. 

b) For <I)o,cr < *^0 < ^max solutious will tend to zero. 

In the absence of ‘Fo.cr only the hrst behaviour holds. For the specihc case $0 = *^ 0,^5 the 
solution will tend to $. So this solution which tends to the saddle point is the separatrix 
which separates the two types of behaviours. Hence for given parameters k, A and c 
satisfying ([T9]), the two types of behaviours are possible provided the existence of the 
separatrix. 

Determination of the initial value $ 0 ,^ yielding the separatrix is a well-known problem 
in the theory of ordinary diherential equations, known as the connexion problem ng. For 
given values of the parameters of the system, it is not possible in general to determine 
analytically whether <Fo,cr exists, neither to hnd a closed form expression for it when it 
does exist. We were able to hnd numerically ‘Fo.cr for certain values of the parameters. For 
instance, for = 20, c = 3 and A = 3, we hnd $o,cr = 0.4936860.... An example of viable 
solutions is shown on Figure |3l 

Using the analytical results given earlier, one can derive an analytical expression for 



( 26 ) 


Let us explain how fl26l) is obtained. As we look for the initial value <Fo,cr for which <F —)■ <h, 
it is clear from flT^ that u{t = 00 ) must then correspond to a zero of the function dn(M | 2 ), 
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Figure 3: The function <F(t) is shown starting at t = 0 at the bounce for the same model 
parameters as in Figure 1. On the left panel, the initial value is <Fo = 0.54 while on 
the right panel $0 = 0.50. As explained in the text <F is positive above $ whose value 
here is <h = 0.40824829..., but can change sign below <F. The curve on the right has a 
more pronounced feature because it starts closer to $o,cr = 0.4936860.... Both curves 
tend asymptotically to zero and constitute viable bouncing solutions. The curve starting 
precisely at <Fo,cr will tend asymptotically to $. 


which are given by 


(2j + 1) F^ (1) 
4^/^ 

Using ffTT|) . we obtain 


CK^ 


<F, 


0 ,cr 


J = 0,1,.... 


Al ^121 + 


Ak2 


x/3 


£!a) 


(27) 


(28) 


Inverting this expression we obtain llOtill 

The value of j must be chosen in such a way that <Fo,cr hes inside the allowed interval 
(IT^ and that fl28|) is satisfied. As the function dn(M|2) is periodic, only a small number 
of values of j must be considered. The numerical value of ^o.cr given above was obtained 
from (126]) for j = 1. For some values of the model parameters and any value of j, ‘Fo.cr is 
not in the allowed interval and then the resulting model is not viable. 

In the presence of radiation one checks easily that the allowed range of initial values 
shrinks with a new larger minimal value <Fq satisfying 

Krmn = ■ ( 29 ) 

We have of course A < A' < 0. The quantities <F, <Fmax on the other hand depend 
on the model parameters and remain unchanged. One can just repeat the proof of the 
existence of viable bouncing solutions. Hence, as long as we have ‘ho^m < ‘^’o.cr, the 
set of viable bouncing solutions remains. In the new expression for <F(f) we will have 
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M dn ^ ^ I 2^ - i\/2 F {i\J^ t \ 2^ while the same corrective 

factor (^)"‘ appears in front of ffl^ as well. We have the obvious corresponding change 
in fl25]) and fl2BD . 

Let us consider next the possibility to have a viable bounce when Z = —1 in the 
absence of radiation. If we want to keep the conformal invariance, the absence of ghosts 
{ubd > —|) and the existence of a bounce, we are led to the following ansatz: 

Fz=_i($) = -Fz=i{^) Uz=-i{^) = -f/z=i($) , (30) 

with the parameters A still positive. The quantities Fz=i, resp. Uz=i, correspond 
to ([5]), resp. ([6]). It is easily checked that ([7]) is retrieved leading to a bounce provided 
the quantity A is negative. This implies here too that c > 0, but now Uz=-i no longer 
represents an inverted potential and we have Uz=-i{^ = 0) < 0. As we have emphasized 
earlier, with (130|) the same equation of motion for the held <I> is obtained therefore leading 
to the same solution (na. However, the crucial difference is that the domain of physical 
validity differs in both cases: now all values d* > ^ satisfy Fz=-i > 0. At the bounce, 
it is seen from (El) that the potential must satisfy Uz=-i{^o) > 0 yielding in this case too 

•ho > ^0,min- 

Up to these differences, we can essentially repeat the analysis performed for Z = 1. It 
is easy to show using (jl]) that $ = 0 implies $ > 0 for $ > $, and $ < 0 for $ < I*. On 
the other hand from ([3]), $ = 0 implies <I> < 0. Hence d) = 0 is not possible for d) > $ 
while it would correspond to a maximum for $ < <h. On the other hand, from (jl]) only 
$ remains an acceptable hnite limit of $ because $ = 0 is now excluded. Imposing the 

A 2 

condition ^ > 1, we now obtain 

36c ’ 

\/6 

- < < $ • (31) 

If <I>o > *h, $0 > 0 implies <I> —)■ cxd while d>o < 0 can in principle produce two possible 
behaviours namely d) —)■ d) (in that case d>o = *ho,cr) or $ ^ 0, the latter being unviable. 
The hrst possibility cannot occur. Indeed with <I)o < 0, some initial values above as well as 
below <ho,cr would tend to zero. But the particular value <ho,cr must separate two different 
behaviours, hence this is not possible. When $o < *^5 0 for a negative initial slope 

because there is no minimum. With <I)o > 0, three different behaviours can be found: 
<h —>■ 0 when $ goes through a maximum (below <h),$—>■ oo when there is no maximum 
and hnally d) ^ $ for <ho,cr which now separates the two behaviours mentioned earlier, its 
existence is therefore allowed in this case. It is interesting that while a viable bounce can 
exist for Z = —1 with a potential bounded from below, it has measure zero in the initial 
conditions, only one particular set of initial conditions at the bounce is acceptable. The 
expression of <ho,cr for this case is also given by fl26|l . 
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4 Conclusions 


We have shown that for Z = 1 only three types of solutions are possible. All stable solutions 
vanishing asymptotically are especially interesting, their initial value $0 at the bounce can 
vary in some range <ho,cr < "ho < The solution with $0 = ‘f’o.cr is unique (measure 

zero) and moreover unstable. The stable solutions lead to the General Relativity limit if 
we assign to the bare gravitational constant its numerical value in GR. Gonsistency of 


our model requires then the inequality 


Stt 


•Cl. If we want to embed our model in a 


realistic cosmology, this attractive feature is even reinforced by the possibility to add a 
radiation component to our bouncing model. 

Further, due to the invariance of the equations under the change t —)■ —t, the solutions 
for f < 0 are obtained from those found for f > 0 changing the sign of <F. Hence, regular 
solutions for t > 0 will diverge in the past both for Z = 1 and Z = —1. Still, starting with 
regular initial values at —00 < f < 0, it is possible to have a bouncing universe regular in 
the future. Note that a solution symmetric around the bounce requires $0 = 0 which will 
diverge both in the past and in the future. 

It is well known that a bounce in the early universe can cure many of the problems 
occurring in Big Bang cosmology besides of course avoiding the Big Bang singularity itself. 
Without exploring further this scenario, we have shown that surprisingly such a bounce 
can occur for a FLRW flat universe within one of the oldest and best understood alter¬ 
natives to GR in a model tending dynamically to GR after the bounce. The key point 
is to admit a negative scalar field potential that can lead to a transient negative energy 
density of the scalar held for some range of parametersjf Of course, to make the model 
cosmologically viable, one has to generalize it somehow to open channels for decay of the 
effective positive (and large) cosmological constant A into other quantum helds including 
those of the Standard Model of elementary particles at late times. 
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